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Abstract 

The ratio of the symmetry energy coefficient to temperature, a sym /T, in Fermi energy heavy ion 
collisions, has been experimentally extracted as a function of the fragment atomic number using 
isoscaling parameters and the variance of the isotope distributions. The extracted values have 
been compared to the results of calculations made with an Antisymmetrized Molecular Dynamics 
(AMD) model employing a statistical decay code to account for deexcitation of excited primary 
fragments. The experimental values are in good agreement with the values calculated but are 
significantly different from those characterizing the yields of the primary AMD fragments. 
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I. INTRODUCTION 



In Fermi energy heavy ion collisions, fragments are copiously produced. The mass distri- 
butions of these fragments exhibit a power low behavior which has been discussed in terms 



of the Modified Fisher Model 



The isotope distributions of these fragments play a key 
role in these analyses. Theoretical studies indicate that the isotope formation is governed 
by the free energy at the density and temperature of the emitting system. The experimental 
observation of isoscaling for two similar reactions with different neutron to proton ratios, 
N/Z, demonstrates that the free energies and therefore the yields of the fragments are also 



closely related to the N/Z of the emitting system 
isotope with N neutrons and Z protons can be given by 
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. Thus the experimental yield of 

mm-. 



Y(N, Z) = Y F(N, Z)A- T exp{-[G(N, Z) 

-fi n N-p p Z]/T} (1) 

where Yq is a constant and G(N,Z) is the nuclear free energy at the time of the fragment 
formation. /i n and \x v are the chemical potential of neutron and proton, and T is the 
temperature of the emitting source. The factor F(N,Z) is the correction factor for the 
feeding from the statistical decay processes. The factor, A~ T , originates from the entropy of 
the fragment jl|. The symmetry energy term in the free energy, G(N,Z), is usually expressed 
as: 

E sym = a sym (N - Zf/A (2) 

where A=N+Z and a sym is the symmetry energy coefficient which depends on the nuclear 
density p and the temperature T of the emitting source. 

From Eq.([T]) R 12 , the ratio of the isotope yields for two similar reaction systems with 
different N/Z ratios, can be written as: 

R 12 (N, Z) = Cexp(aN + fiZ) (3) 

This relation is known as the isoscaling relation. The isoscaling parameters, a = {p\— /i^) /T 
and j3 = {l^—pfy/T) are the differences of the neutron or proton chemical potentials between 
the systems 1 and 2, divided by the temperature. C is a constant. System 1 is normally 
taken as the more neutron rich of the two. 
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As discussed in refs. 0, ^ [9], the isoscaling parameters and the symmetry energy coef- 
ficient are closely related. For a multifragmentation regime, as pointed out in ref. [9(, this 
relation is given by: 

a(Z) = Aa sym A(Z/A) 2 /T (4) 

where A(Z/A) 2 = (Z/A) 2 1 — (Z/A) 2 2 for the two reaction systems and A is the average 
mass number of isotopes for a given Z. 

There are two issues for the determination of the a sym values in Eq.(j3J). One is the 
source temperature T. Since the beginning of the experimental study of heavy ion collisions 
in the multifragmentation regime, significant efforts have been made to extract the source 
temperature, but different methods of temperature extraction can lead to different results 
and uncertainties still remain [h]]. Another issue is the effect of the secondary decay process, 
expressed by F(N,Z) in Eq.flTJ). In experiments fragments have typically cooled down to the 
ground state before they are detected. Indeed, in previous works, excitation energies of the 
primary fragments have been evaluated by studying the associated light charged particle 
multiplicities Q. Such data raise the question of the degree of confidence for the 
experimentally extracted symmetry energy coefficient, in which this important effect is not 
properly corrected. In fact in a separate paper using the same data set presented here we 
have demonstrated that the secondary decay processes significantly effect the isobaric yield 
ratios and the experimentally extracted symmetry energy coefficient [13]. In this paper, we 
focus on the relation between the ratio a sym /T and isoscaling parameters and between the 
ratio and the widths of the isotope distributions. The experimentally extracted values of 
a sym /T extracted from both observables are compared to those extracted from the model 
calculations. 



II. EXPERIMENT 



The experiment was performed at the K-500 superconducting cyclotron facility at Texas 
A&M University. 64,70 Zn and 64 Ni beams were used to irradiate 58 ' 64 Ni, 112 ' 124 Sn, 197 Au and 
232 Th targets at 40 A MeV. Intermediate mass fragments (IMFs) were detected by a detector 
telescope placed at 20°. The telescope consisted of four Si detectors. Each Si detector was 
5cm x 5cm. The nominal thicknesses were 129, 300, 1000, 1000 fim. All Si detectors 
were segmented into four sections and each quadrant had a 5° opening angle in polar and 
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azimuthal angles. Therefore the energies of the fragments were measured at two polar 
angles of the quadrant detector, namely 9 = 17.5° ± 2.5° and 9 = 22.5° ± 2.5°. Typically 
6-8 isotopes for a given atomic numbe up to Z=18 were clearly identified with the energy 
threshold of 4-10 A MeV, using the AE-E technique for any two consecutive detectors. The 
AE-E spectrum was linearized empirically. Mass identification of the isotopes were made 
using a range-energy table [15[. In the analysis code, isotopes are identified by a parameter 
ZreoI- For the isotope with A=2Z, ZR ea i = Z is assigned and other isotopes are identified by 
interpolation between them. Typical Z^ ea i spectra are shown in Fig.l. The energy spectrum 
of each isotope was extracted by gating the isotope in a 2D plot of Zn ea i vs energy. The 
yields of light charged particles (LCPs) in coincidence with IMFs were also measured using 
16 single crystal CsI(Tl) detectors of 3cm thickness set around the target. The light output 
from each detector was read by a photo multiplier tube. The pulse shape discrimination 
method was used to identify p, d, t, h and a particles. The energy calibration for these 
particles were performed using Si detectors (50 -300 /^m) in front of the Csl detectors in 
separate runs. 

The yield of each isotope was evaluated, using a moving source fit. For LCPs, three 
sources (projectile-like(PLF), nucleon-nucleon-like(NN) and target-like (TLF)) were used. 
The NN-like sources have source velocities of about a half of the beam velocity. The pa- 
rameters are searched globally for all 16 angles. For IMFs, since the energy spectra were 
measured only at the two angles of the quadrant detector, the spectra were parameterized 
using a single NN-source. Using a source with a smeared source velocity around half the 
beam velocity, the fitting parameters were first determined from the spectrum summed over 
all isotopes for a given Z, assuming A=2Z. Then all extracted parameters except for the nor- 
malizing yield parameter were used for the individual isotopes. This procedure was based 
on the assumption that, when the spectrum is plotted in energy per nucleon, the shape of 
the energy spectrum is the same for all isotopes for a given Z. Indeed the observed energy 
spectra of isotopes are well reproduced by this method. For IMFs, a further correction was 
made for the background. As seen in Fig.l, the isotopes away from the stability line, such 
as 10 C and 29 Mg, have a very small yields and the background contribution is significant. 
In order to evaluate the background contribution to the extracted yield from the source fit, 
a two Gaussian fit to each isotope combined with a linear background was used. The fits 
are shown in Fig.l. Each peak consists of two Gaussians. The second Gaussian (about 10% 
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of the height of the first one) is added to reproduce the shape of the valley between two 
isotopes. This component is attributed to the reactions of the isotope in the Si detector. 
The centroid of the Gaussians was set to the value calculated from the range-energy table 
within a small margin. The final yield of an isotope was determined by correcting the yield 
evaluated from the moving source fit by the ratio between the two Gaussian yields and the 
linear background. Rather large errors ( ~ ±10%) are assigned for the multiplicity of the 
NN source for IMFs, originating from the source fit besides the background estimation. The 
errors from the source fit are evaluated from the different assumptions of the parameter set 
for the source velocity and temperature. 
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FIG. 1: Typical linearized isotope spectra are shown for Z=6 (upper figure) and 12 (lower) cases. 
The number at the top of each peak is the mass number assigned. Linear back ground is assumed 
from valley to valley for a given Z. Each Gaussian indicates the yield of the isotope above the back 
ground. 
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III. ISOSCALING 



At the top of Fig. [21 the yield ratio of Eq.© for the reactions 64 Ni + 124 Sn and 64 Ni + 58 Ni 
is plotted as a function of N. The a parameter is determined by individual fits to yield ratios 
for isotopes with a given Z. The extracted values are plotted in the bottom of Fig. [2J As seen 
in the figure, the extracted ot(Z) parameter shows a steady decrease as Z increases for Z > 
4. The (3{N) parameter generally shows much similar variation with increasing N, and has 
the opposite sign. Hereafter a and are denoted as a(Z) and 0(N). Isoscaling parameters, 
oi{Z) and (3(Z), have been evaluated for all possible combinations between two reactions. 
For 18 different reactions are considered here. More than 150 different combinations have 
been studied. 
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FIG. 2: (upper) a(Z) values as a function of Z for 64 Ni + 124 Sn and 64 Ni + 58 Ni. From the left to 
right each lines correspond to Z=l to Z=18. (lower) The extracted a(Z) values are plotted as a 
function of Z. 



In Fig. [3] the extracted a>(Z) values are plotted as a function of A(Z/A) 2 for Z=6 and 
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Z=12. Each data point represents a combination of two reactions. As seen in the figure, 
the oc(Z) values are linearly related to A(Z/A) 2 . The slope associated with this relationship 
increases gradually as Z increases. The correlations have been fit by a linear function for 
each Z and the slope values, which correspond to the value 4a sym /T in Eq.(T4]), have been 
extracted. 

In Fig. [4] the extracted values of a sym /T are plotted as a function of Z and shown by solid 
circles. A clear trend is observed for the parameter, a sym /T. The value increases from 4 to 
14 as Z increases from 4 to 15. The value for Z=3 is much larger than Z=4. This is partially 
caused by the isotope distribution. Since 5 Li is unstable and decays before arriving to the 
detector, A deviates from the actual centroid of the isotope distribution. An attempt has 
been made to fit the distribution by a Gaussian function and determine A as the centroid 
value. This procedure makes a sym /T around 3 for Z=3, but the uncertainty is significant, 
especially for neutron deficient systems. Therefore in the plot the experimental Z/A value 
appears without correction. 
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FIG. 3: a(Z) values as a function of A(Z/A) 2 for Z=6 (upper) and Z=12 (lower). The lines are 
the results of a linear fit. 
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FIG. 4: Experimental a sym (p,T)/T values extracted from Eq.Q are shown by solid circles as a 
function of Z. Open circles show calculated a sym (p,T)/T values for the fragments yields in AMD- 
Gemini calculations. Squares are experimental results of Q{Z) discussed in the section IV. 

IV. SYMMETRY ENERGY AND VARIANCE OF THE ISOTOPE DISTRIBU- 
TION 

The multiplicity distributions of the isotopes for a given Z show a quadratic distribution 
when they are plotted on a logarithmic scale. Since the symmetry energy term is the 
only term proportional to (N-Z) 2 in the free energy, this suggests that the variance of the 
distributions is closely related to the symmetry energy coefficient. In order to explore the 
relation between the symmetry energy term in the free energy and the variance of the 
isotope distribution, Ono et al. introduced a generalized function K(N,Z) for the free energy 
in ref. Q as given below. 

n 

K(N, Z) = J2 u>i(N, Z)[-lnYi(N, Z) + a t (Z)N + 7i (Z)] (5) 
i=i 

Here i represents each reaction. The summation is taken over i for the different N/Z reaction 
systems in order to get isotope multiplicity distribution in a wide range from proton rich 
to neutron rich isotopes. The average weights, Wj(N,Z), are determined by minimizing 
the statistical errors in K(N,Z) for a given (N,Z). The isoscaling parameter, a^Z\ is the 
isoscaling parameter value evaluated in the previous section. For each Z the parameters, 
7i(Z), are determined by optimizing the agreement of the quantities [-lnY(N,Z) + a.i(Z) + 
7i(Z)] from different reactions. A typical K(N,Z) distribution from the experiment is shown 
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in Fig. [5j The isotope distributions for a given Z exhibit a smooth quadratic distribution 
and they can be fit by a function: 

K(N, Z) = £(Z)N + V (Z) + C(Z)(N - Zf/A (6) 

Where £,{Z), 1](Z), ((Z) are the fitting parameters. As one can see the functional form, £(Z) 
is related to the symmetry energy coefficient given in Eq.(j3|) as: 

((Z) = a sym /T (7) 

In Fig. HI the values of £(Z) extracted using this technique are shown by solid squares. The 
values are generally about 1 or 2 units smaller than the a sym /T values evaluated in the 
previous section (solid circles), but the general trend is in good agreement. The difference in 
the extracted values for Z=3 originate from the different determination of the average value 
of A. In the analysis in this section the average masses of these isotopes are determined by 
the centroid of the quadratic distributions. 




V. COMPARISONS WITH MODEL SIMULATIONS 

The experimentally detected fragments are the final products of the reaction. Excited 
primary fragments will have cooled down by statistical decay before they arrive in the 
detector. Excitation energies of the primary fragments have been evaluated experimentally 
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by measuring the light charged particle multiplicities in coincidence with the fragments 



121 ] . Typical excitation energies of 2 to 3 MeV/nucleon have been derived. In order to study 



the effect of the secondary decay process on the experimentally extracted ratio, a. ^ /T , the 
simulation codes of an Antisymmetrized Molecular Dynamics (AMD) model [l6|, ll7|, |24| and 
a statistical decay code, Gemini [19|, have been used. These codes have often been used to 
study the fragment production in Fermi energy heavy ion reactions and the global features of 
the experimental results have been well reproduced jl^. I20I 24 ] . Since the AMD calculation 
requires a lot of CPU time, only two of the experimental reaction systems have been studied. 
The systems examined are 64 Zn+ 112 Sn and 64 Ni+ 124 Sn at 40 AMeV. All calculations shown 
in this pa per have been performed in a newly installed computer cluster in the Cyclotron 
Institute [25]. In order to obtain yields of the final products, the deexcitation of primary 
fragments formed at 300 fm/c was followed using the Gemini code until they cooled to the 
ground state. Using the same analysis described in the section III, the scaling parameters, 
ce(Z) and /3(Z) were then extracted from the simulated events as a function of Z. The average 
mass number of the isotopes for a given Z was also evaluated from the calculations. The 
calculated variation of a sym /T is shown by open circles in Fig. HI The values are typically 
one to two units higher than the experimental values (open circles) but exhibit an essentially 
identical trend. For the primary fragments at t = 300 fm/c the same analysis has been made. 
These fragments are identified using a coalescence technique in phase space. The evaluated 
symmetry coefficients using R c = 5 are shown in Fig. [6l R c =5 corresponds to a radius of 
5 fm in configuration space. For the primary fragments the calculated values are close to 
the experimental values observed for Z <6 but, in contrast to the experimental results, they 
remain more or less constant for all Z. The rather flat distribution of a sym /T values over the 
entire range of Z is consistent with the picture of the origin of the primary fragments from a 
common emitting source with a given density and temperature. These comparisons indicate 
that the a sym /T values are significantly modified by the secondary process for Z > 4 



VI. SUMMARY 



The symmetry energy coefficient to temperature ratio, a sym /T , as a function of Z has 
been extracted from experimental data in two different ways, namely from isoscaling pa- 
rameters and from the variance of the observed isotope distributions. The results from the 
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FIG. 6: Comparisons between the experimentally extracted symmetry energy coefficient to tem- 
perature ratio and those from AMD. Solid circles show the experimental values determined from 
Eq.(H]) and open circles show results from of AMD-Gemini calculations. Closed squares indicate 
results for the primary fragments with Rc=5. 

two techniques are in reasonable agreement. Experimental values increase from ~ 4 to 12 
as Z increases from 4 to 15. The values and trends observed for the final fragments are 
well reproduced by the AMD plus Gemini model simulations. However these values are 
significantly different from the values extracted for the primary fragments, especially for Z 
> 4, suggesting that the derivation of the observed ratio of symmetry energy coefficient to 
temperature is significantly perturbed by the secondary decay processes. 
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